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“Unsteady  and  Separated  Flows” 
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Research  on  "Unsteady  and  Separated  Flows” 


The  ultimate  aim  of  our  research  is  to  obtain  understanding  of  and  develop  methods  for  analyzing - 

and  for  controlling  unsteady  separated  flows;  included  are  flows  which  are  globally  steady  but  have - ^ 

unsteadiness(usually  periodic)  due  to  instability  of  the  separated  flow,  as  well  as  flows  which  are  intrin- 


sicalty  unsteady  due  to  unsteady  boundary  conditions,  either  periodic  or  impulsive.  Our  approach  is  by 
experiments  in  wind  tunnels  and  water  facilities  as  well  as  by  numerical  simulation.  A  large  part  of  the 
effort  includes  the  acquisition  and  development  of  tools  and  techniques  for  the  laboratory  experiments 
and  the  development  of  computational  methods. 


□ 


The  following  report  is  divided  into  three  parts,  based  on  the  areas  of  emphasis  of  the  three  principal  Cod#8 

lAvo?!.  and/or 


investigators;  the  activities  in  these  areas  are  in  fact  overlapping  and  interactive. 
A.  Globally  Steady  and  Intrinsically  Unsteady  Separated  Flows 
1.  Research  Objectives 


Dist 


m 


Special 


The  objective  of  this  part  of  the  research  is  to  investigate  the  evolution  of  separated  flow  in 
unsteady,  mainly  impulsive  motions  of  simple  but  airfoil-related  bodies.  This  evolution  involves  the 
shedding  of  vorticity,  the  development  of  the  vorticity  into  large  vortical  structures  and  their  interaction 
with  each  other  and  with  the  body.  An  important  task  is  to  to  correlate  the  flow  development  with  the 
history  of  forces  (lift,  drag  and  moment)  on  the  body.  Accordingly  our  development  of  experimental 
facilities  and  techniques  has  been  aimed  at  flow  field  observation  and  measurement  as  well  as  measure¬ 
ment  of  forces. 


2.  Status  of  Research 

\ 

(a)  Investigations  in  the  X-Y  Tow  Tank 

While  some  unsteady,  non  periodic  maneuvers,  such  as  pitch  up,  can  be  investigated  in  steady-flow 
facilities  (wind  or  water  tunnels),  unsteady  trajectories  in  general  require  a  special  facility  such  as  our 
X-Y  Tow  Tank.  A  large  part  of  this  year's  effort  has  been  the  continuing  development  of  experimental 
techniques  for  this  facility.  Figures  1  and  2  and  the  following  paragraphs  summarize  these  develop¬ 
ments. 


X-Y  Carriage  Computer  Control; 


The  computer  control  system  for  the  X-Y  carriage  was  completed  and  installed  or,  the  ne*v  XT  com¬ 
putable  microcomputers.  In  this  system,  a  GalicI  position  control  board  receives  position  feedback  from 
optical  encoders  on  the  carriage  drive  motors  and  generates  amplifier  command  voltages.  Interactive 
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software  is  used  to  generate  real  time  incremental  position  commands  for  the  Galiel  board;  the  program 
was  designed  to  generate  complex  (unsteady)  carriage  motion  profiles  in  a  user  friendly  manner.  Some 
final  adjustment  of  the  control  system  parameters  is  still  required  to  optimize  carriage  response  vs. 
bibration  level.  This  is  of  immediate  concern  for  force  measurements  since  the  very  slight  vibration 
induced  by  the  position  control  system  is  apparent  in  the  force  output  traces. 

Load  Cell  Force  Balance; 

This  system  consists  of  three  components:  the  force  balance  itself,  which  is  designed  to  be 
suspended  below  the  X-Y  carriage  and  will  measure  3  force  components  (lift,  drag  and  moment  in  a 
horizontal  plane,  a  four  channel  amplifier/analog  noise  filter;  and  an  RC-electronics  computerscope  data 
acquisition  system  running  on  an  AT  compatible  computer.  Output  data  are  analyzed  and  optimally 
filtered  using  software  which  can  output  either  forces  or  force  coefficients  vs.  time  and/or  position. 
Amplifier  gain  can  be  adjusted  to  allow  full  scale  forces  of  ±  20N  to  ±  .02N,  the  amplifier  noise  being 
about  .002N.  A  4th  order  low  pass  25  Hz.  filter  is  used  to  analog  filter  the  signal  before  it  is  digitized 
to  eliminate  aliasing  in  any  subsequent  digital  or  Fourier  filtering. 

The  balance  was  designed  to  have  a  minimum  natural  frequency  of  25  Hz.,  but  initial  force 
measurements  indicated  the  very  slight  vibration  of  the  carriage  due  to  the  computer  control  system 
induced  a  vibration  somewhere  in  the  balance  system  of  about  8  Hz.  This  resulted  in  a  signal  to  noise 
ratio  of  about  2.  A  little  further  work  should  eliminate  this  vibration  and  a  signal  to  noise  of  30+ 
should  be  achievable. 

Preliminary  force  measurements  have  been  made  on  a  specially  constructed  flat  plate  set  at  different 
angles  of  attack,  using  different  acceleration  profiles  (both  linear  and  sinusoidal)  and  at  different 
Reynolds  numbers.  Zero-acceleration  profiles  (i.e.,  globally  steady  motion)  were  also  investigated.  In 
the  latter  case,  with  the  flat  plate  at  angle  of  attack  90°,  drag  histories  such  as  that  shown  in  Figure  3 
were  obtained.  Although  the  nearly  steady-state  part  of  the  history  is  at  Cd  =  2,  close  to  the  nominal 
value  for  flat-plate  drag,  there  are  some  puzzling  features,  for  example  the  early  onset  (x/c  =  3)  of  this 
"steady"  state,  as  well  as  the  not  quite  steadiness  (before  termination  at  x/c  =23).  Comparison  with 
our  numerical  computations  are  discussed  in  the  following  section  C. 

(b)  Investigations  in  Wind  and  Water  Tunnels 

The  earlier  steady  state  experiments  in  tow  tank  and  wind  tunnel  on  a  flat  plate  with  vortex 
"trapped"  ahead  of  an  upper  flap  (the  Hurlcy-Tanveer-Saffinan  configuration)  (Figure  4  and  5)  showed 
that  the  full  theoretical  lift  cannot  be  achieved  because  the  flow  separates  from  the  upper  flap  soon  after 
reattaching  to  its  leading  edge  (due  to  the  adverse  pressure  gradient  on  the  upper  surface  of  the  upper 
flap).  Various  combinations  of  plate  lengths,  angle  between  plates,  and  angle  of  attack  were 
investigated,  as  well  as  addition  of  curvature  to  the  uper  plate.  Some  combinations  of  parameters 
showed  lift  coefficients  as  high  as  2.2  at  large  angles  of  attack  (a  =  40°),  even  with  large  separation.  It 
was  thought  that  by  properly  designing  the  upper  plate,  using  techniques  suggested  by  Stratford,  the 
recovery  on  the  upper  plate  could  be  enhanced.  This  led  to  the  idea  of  using  a  Liebeck  airfoil  (which  is 
designed  to  recover  upper  surface  pressure)  with  a  large  cutout  on  that  part  of  the  upper  surface  which 
is  at  constant  Cp  (corresponding  to  the  free  streamline  in  the  Hurley-Tanveer-Saffman  model).  We 
hoped  that  the  cavity’s  free-shear-layer  oscillation  would  provide  the  strong  vortices  necessary  to 
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energize  the  boundary  layer  ahead  of  the  adverse  pressure  gradient  region,  as  a  past  experiment  by 
Gharib  has  shown. 

We  selected  and  constructed  a  Liebeck  airfoil  that  shows  the  highest  theoretical  lift  coefficient 
available  (CL  -  8)  in  the  Liebeck  airfoil  series.  Tests  were  conducted  on  a  large  (1.5  ft  chord  length 
and  4  ft  span)  model  in  a  water  tunnel  at  the  University  of  Southern  California,  San  Diego.  Initial  flow 
visualization  studies  indicate  that  the  presence  of  the  cavity  improves  the  reattachment  length  (Figure 
5).  The  extent  of  the  reattachment  region  is  a  strong  function  of  flow  condition  in  the  cavity.  A  variety 
of  cavity  widths  were  tested  for  the  optimum  and  non-optimum  angles  of  attack.  It  was  concluded  that 
the  condition  of  self-sustained  oscillations  in  the  cavity  is  very  effective  in  extending  the  reattachment 
region,  but  depends  strongly  on  the  angle  of  attack.  The  optimum  angle  of  attack  for  the  selected 
Liebeck  airfoil  has  a  narrow  range.  As  a  consequence,  a  new  Liebeck  airfoil  (LA  505  series),  with 
milder  adverse  pressure  gradient  and  broader  range  of  angle  of  attack  for  optimum  lift,  has  been 
designed.  It  also  became  clear  that  measurement  of  lift  and  drag  is  important  for  evluating  such 
configurations;  consequently  a  force  balance  for  these  experiments  has  been  designed  and  will  be  built 
(with  UCSD  funds).  The  Liebeck  airfoils  and  the  force  balance  are  designed  to  be  compatible  with  both 
the  UCSD  water  tunnel  and  the  GALCTT  Merrill  Wind  Tunnel. 

Work  is  still  needed  to  improve  the  quality  and  quantity  of  the  output  traces. 

Particle  Tracing  (Streakline)  Flow  Visualization  System: 

The  particle  tracing  system  consists  of  an  American  Laser  Corporation  8W  Argon-Ion  laser  whose 
intensity  is  controlled  using  a  Bragg  cell.  The  laser  beam  is  split  into  a  sheet  and  is  used  to  illuminate 
fluorescent  tracer  particles  in  the  flow;  output  traces  are  then  recorded  with  an  integrating  CID  camera. 
The  time  of  intergration  and  the  pulsation  of  the  laser  beam  through  the  Bragg  cell  is  controlled  using  a 
two  channel  programmable  function  generator  which  is  controlled  by  an  XT  compatible  computer. 
Final  video  images  are  digitized  using  a  PCVision  image  capture  board  operating  in  an  AT  compatible 
computer,  analysis  software  is  then  used  to  enhance  the  video  image  and  calculate  velocity  and  vorticity 
fields. 

Work  is  still  needed  to  improve  the  quality  and  quantity  of  the  output  traces.  The  system  in  its 
present  form  uses  a  cylindrical  lens  to  create  a  sheet  of  laser  light.  Since  the  trace  particles  fluoresce 
more  efficiently  at  higher  incident  light  levels,  using  a  scanning  mirror  to  create  the  sheet  should 
improve  the  sharpness  of  the  particle  streaks;  the  particles  would  then  be  "hit"  with  the  full  laser  power 
in  pulses,  instead  of  a  lower  level  steady  lighting.  Some  development  work  will  be  required  for  the 
technique  to  be  useful  on  unsteady  or  "complex"  flows  (i.e.,  a  flow  field  with  many  small  structures 
present)  since  most  of  its  previous  development  has  been  with  steady  flows  where  the  magnitude  of  the 
velocity  is  relatively  constant  and  the  particle  density  can  be  fairly  low. 

B.  Feedback  Control  of  Separated  Flows 

1.  Research  Objectives 

The  objective  of  this  research  is  to  apply  feedback  control  to  separated  shear  flows.  Due  to  the  very 
high  order  and  nonlinear  nature  of  this  control  problem,  our  approach  will  be  to  use  an  "adaptive” 
control  method.  We  plan  to  integrate  simplified  equations  of  motion  for  the  flow  system  in  real  time. 
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These  equations  can  be  thought  of  as  the  system  "model".  Parameters  for  this  system  will  be  adjusted 
in  response  to  control  inputs,  real-time  flow  measurements,  and  by  direct  integration.  For  example,  if  a 
vortex-dynamics  calculation  were  used,  the  parameters  would  be  the  locations  and  strengths  of  the 
vortices.  The  feedback  loop  will  be  closed  by  continuously  trying  to  decrease  a  function  such  as  the 
energy  in  the  wake  (a  candidate  Liapunov  function)  with  respect  to  possible  control  inputs.  Initially, 
the  combined  flow  and  controller  will  be  simulated  without  the  actual  flow,  to  determine  the 
effectiveness  of  this  approach. 

2.  Status  of  Research 

Exploratory  open  loop  experiments  have  been  performed  on  smooth  circular  cylinders  in  steady 
uniform  flow,  executing  forced  rotary  oscillations.  Our  flow  visualization  and  wake  profile 
measurements  at  Reynolds  numbers  in  the  vicinity  of  16,000  have  shown  that  it  is  possible  to  produce 
dramatic  effects  on  the  structure  of  the  cylinder  wake  by  such  means.  In  particular,  a  factor  of  five 
decrease  in  the  resulting  displacement  thickness  and  estimated  cylinder  drag  was  achieved,  depending 
on  the  frequency  and  amplitude  of  oscillation. 

The  cylinder  was  rotated  sinusoidally  at  a  forced  Strouhal  number  (S/)  of  fd/U frte_strtam,  where 
/  =  forcing  frequency,  d  =  cylinder  diameter,  and  Ufret-siream  =  free  stream  velocity,  while  the  peak 
circumferential  velocity  (vmax)  was  fixed  at  2 Ufrte-Slream.  This  value  of  vm„  was  chosen  because  it 
is  greater  than  the  maximum  flow  velocity  in  the  vicinity  of  the  unforced  cylinder.  As  such,  we 
expected  that  at  the  peak  of  the  oscillation  cycle,  only  one  sign  of  vorticity  would  be  created  on  the 
cylinder  surface.  The  wake  mean  velocity  and  rms  profiles  reported  in  Figure  6  were  measured  four 
diameters  downstream  of  the  cylinder.  Of  particular  interest  is  the  coincidence  of  observations  in  the 
neighborhood  of  Sf  -  0.9.  Here  we  found  minima  for  the  drag  coefficient  (CD)  and  the  displacement 
thickness  (8/d),  a  transition  from  forcing  of  the  wake  vortices  to  forcing  of  shear  layer  vortices,  and  the 
ratio  of  vertical  to  horizontal  vortex  spacing  (a)  was  of  order  one,  transitioning  from  a  <  1  to  c  >  1 
with  increasing  5/ . 

It  was  observed  over  the  entire  range  of  observation  (0  <Sf  <  2.5)  that  the  vortex  shedding  was 
synchronized  with  the  forcing  frequency.  Visually,  the  wake  structure  took  on  three  basic  personalities. 
When  Sf  «.  0.2  it  was  observed  that  the  flow  in  the  vicinity  of  the  cylinder  can  look  somewhat  like 
the  potential  flow  around  a  cylinder  with  circulation.  In  this  region,  the  flow  was  characterized  by  the 
cylinder  swinging  the  flow  back  and  forth  via  the  Magnus  effect.  Large  vortices  of  alternating  sign 
were  shed  between  the  positive  and  negative  parts  of  the  forcing  cycle.  As  Sf  increases  past 
Sf  =  0.2,  the  method  of  flow  actuation  changed,  vorticity  was  collected  on  the  low-speed  side  of  the 
cylinder,  while  being  released  on  the  high-speed  side.  Here  the  "Magnus  effect"  was  not  apparent,  and 
the  flow  around  the  cylinder  was  characterized  by  the  shed  vortices.  The  resulting  patterns  of 
alternating  vortices  persisted  far  downstream  of  the  cylinder  for  Sf  <  0.8.  In  contrast,  for  Sf  >  0.9, 
only  the  shear  layers  separating  from  the  cylinder  were  forced;  these  two  rows  of  shear  layer  vortices 
were  drawn  toward  each  other,  nearly  closing  the  wake.  They  disintegrated  upon  meeting,  with  the 
wake  transitioning  at  a  lower  frequency  soon  afterward. 

The  displacement  thickness  and  drag  coefficient  estimates  summarized  in  Figure  7  were  calculated 
from  the  mean  velocity  profiles  in  Figure  6.  Note  the  broad  minima  in  Cp  and  5 Id  around  Sf  =  0.9. 
There  may  also  be  one  or  more  local  minima  below  Sf  =  0.3,  but  this  possibility  was  not  investigated. 
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The  validity  of  the  estimated  CD  below  Sf  -  0.2  is  in  question  because  of  the  large  velocity 
fluctuations.  We  suspect  that  Cq  has  been  underestimated  in  that  region.  The  error  bars  in  Figure  7b 
reflect  the  uncertainty  due  to  the  velocity  fluctuations. 

These  open  loop  experiments  suggest  that  control  of  bluff  body  flows  can  be  accomplished  by 
encouraging  the  inevitable  (alternating)  vortices  to  occur  at  a  specified  strength  and  rate,  using  the 
subsequent  vortex  interactions  to  control  the  wake. 

C.  Numerical  Simulation  of  Separated  Flows 

1.  Research  Objectives 

This  work  has  aimed  to  perform,  using  vortex  methods,  two-  and  three-dimensional  numerical 
simulations  of  the  dynamics  of  large-scale  separated  vortices  in  the  near  wake  of  a  bluff  body  and  to 
study  the  mutual  interactions  between  the  vortex  structures  and  the  body  during  transient  motions.  A 
primary  objective  in  this  part  of  the  research  is  to  complement  and  interact  strongly  with  the 
experimental  parts  of  the  research.  Our  objectives  for  the  development  of  the  required  computational 
tools  suitable  computational  elements  and  dynamical  equations,  include  the  following  capabilities: 


(a)  Represent  the  vortex  dynamics  in  the  near  and  far  wake  using  suitable  computational 
elements  and  dynamical  equations. 

(b)  Satisfy  the  inviscid  boundary  condition  at  the  solid  surface,  i.e.,  no  flow  through  the  walls, 
and 

(c)  Determine  the  locations  of  the  separation  points  in  time,  the  flux  of  vorticity  into  the  outer 
flow  and,  if  necessary,  faithfully  represent  the  state  of  the  boundary  layer  in  the  vicinity  of 
reattachment. 

2.  Status  of  Research 

(a)  To  represent  the  vortex  dynamics  in  the  near  and  far  wake  of  a  three-dimensional  separated  flow,  we 
have  developed  a  robust  three-dimensional  vortex  method  using  vector-valued  vortex  particles.  These 
vectors  move  with  the  local  velocity  and  stretch  and  rotate  in  response  to  the  local  strain  rate.  A  major 
problem  that  was  encountered  was  the  treatment  of  dissipation.  Without  some  form  of  dissipation  or 
subgrid  scale  modeling,  vortex  calculations,  which  attempt  to  satisfy  the  inviscid  equations  of  motion, 
typically  produce  singularities  in  finite  time  that  terminate  the  simulation.  We  have  derived  a  nonlinear 
viscous  model  for  vortex  simulations  that  is  related  to  the  Smagorinsky  model  for  grid-based  large-eddy 
simulations.  The  method  was  tested  successfully  on  two  separate  problems  involving  the  interaction  of 
two  vortex  rings.  Several  notable  phenomena  were  observed  during  these  interactions:  (1)  complex 
motions  of  the  space  curves  of  die  vortex  tubes,  (2)  development  of  internal  structure  (other  than 
swirling  flow)  such  as  axial  flows  and  vortex  breakdown,  and  (3)  reconnection  of  vortex  lines. 

(b)  In  the  numerical  simulation  of  unsteady  and  separated  flows  of  interest  we  use  vortex  panels  to 
represent  the  boundary-layer  vorticity.  The  usual  formulation  for  the  panel  method  on  a  flat  plate  leads 
to  a  singularity.  Recently,  we  have  derived  a  new  algorithm  to  determine  the  local  boundary  layer 


circulation  per  unit  area  on  either  side  of  a  flat  plate  given  an  arbitrary  wake  vorticity  distribution.  This 
allows  us  to  construct  boundary-layer  separation  models  that  have  a  firm  physical  basis.  We  are  now 
testing  the  new  scheme  on  a  variety  of  two-dimensional  flows  past  a  flat  plate  including  oscillations  in 
angle  of  attack  with  very  encouraging  results.  In  some  cases  (e.g.,  drag  of  a  flat  plate  at  90°  angle  of 
attack)  the  aerodynamic  forces  are  higher  than  those  measured  in  various  experiments  including  the 
present  ones.  Clearly  more  understanding  and  refinement  of  the  experiments  and  the  computations  are 
needed.  We  are  planning  further  computations  and  tow-tank  experiments  to  determine  if  the 
discrepancy  is  due  to  a  modeling  or  numerical  problem  in  the  computation  or  three-dimensional  effects 
in  the  experiments. 

(c)  As  a  first  application  of  our  computational  tools  developed  for  two-dimensional  flows,  we 
investigated  flow  past  a  flat  plate  for  a  fixed  90°  angle  of  attack  and  for  sinusoidal  oscillations  in  angle 
of  attack  about  90°.  For  a  flat  plate  held  at  a  constant  a  =  90°  we  obtained  the  time-dependent  drag 
coefficient  ( CD )  shown  in  Figure  8.  During  the  first  25  units  of  time,  i.e.,  xlc  =25  ({/„.  =  1  and 
plate  length  L  =  1)  the  flow  is  unsteady  but  nearly  symmetric  about  the  horizontal  axis  (see  Figure  9) 
and  CD  drops  below  1.0,  consistent  with  free-streamline  theory.  However,  for  t  >  30,  the  flow 
becomes  highly  asymmetric  and  evolves  into  the  expected,  nearly  periodic  vortex  shedding  mode.  The 
average  drag  is  high  ( <CD>  -  3.6)  as  a  strong  vortex  is  formed  near  each  end  of  the  plate  (Figure  10) 
and  remains  there  for  a  substantial  fraction  of  each  half-cycle,  resulting  in  a  high  average  suction 
distribution  on  the  backside  of  the  plate. 

Inspired  by  cylinder  wake  control  experiments  performed  in  the  experimental  part  of  the  research,  in 
which  rotary  oscillations  of  the  cylinder  dramatically  reduced  the  drag,  we  investigated  the  effect  of 
angle  of  attack  oscillations  on  the  drag  of  a  flat  plate  with  an  average  angle  of  attack  of  90°.  Using  an 
oscillation  frequency  of  about  five  times  the  natural  shedding  frequency  and  a  maximum  tip  speed  of 
2.0,  we  found  a  significant  drop  in  the  average  drag  to  <Co>  -  2.3  (Figure  11)  along  with  significant 
change  in  the  flow  structure.  Under  the  new  conditions  with  oscillations,  two  vortices  are  formed  at 
each  end  of  the  plate  during  each  half  cycle.  The  two  vortices  have  opposite  sign  circulation  and 
therefore  induce  each  other  to  move  away  from  the  plate  (Figure  12)  so  that  the  suction  on  the  backside 
of  the  plate  is  reduced. 
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A.  Roshko  was  Senior  Visiting  Fellow  at  the  Stanford/Ames  Center  for  Turbulence  Research  during 
April  and  May  of  this  year. 

Capt  John  Wissler  of  the  U.S.  Air  Force  Academy  is  enrolled  as  a  graduate  student  in  Aeronautics 
here.  He  will  soon  be  deciding  on  a  research  problem  in  some  area  of  unsteady  flow  and  its  control. 
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Discussions  and  other  communication  h$  continued  with  Dr.  M.  Koochesfahani,  presently  Assistant 
Professor,  Michigan  State  University. 


